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Sufficient conditions for complete controllability of N-level quantum systems subject to a single 
control pulse that addresses multiple allowed transitions concurrently are established. The results 
are applied in particular to Morse and harmonic oscillator systems, as well as some systems with 
degenerate energy levels. Controllability of these model systems is of special interest since they have 
many applications in physics, e.g., Morse and harmonic oscillators serve as models for molecular 
bonds, and the standard control approach of using a sequence of frequency-selective pulses to address 
a single transition at a time is either not applicable or only of limited utility for such systems. 
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I. INTRODUCTION 

Recent advances in laser technology have opened up 
new possibilities for laser control of quantum phenom- 
ena such as control of molecular quantum states, chem- 
ical reaction dynamics or quantum computers. This 
has prompted researchers to study these systems from a 
control-theoretical point of view, in particular in view of 
the limited success of initially advocated control schemes 
based largely on physical intuition in both theory and 
experiment [Xj. 

One issue that arises is the question whether, or un- 
der which conditions, it is possible to control a quan- 
tum system in such a way as to achieve any physically 
permitted evolution of the system. Complete controlla- 
bility is an important theoretical concept that also has 
significant practical implications. For example, it has 
been shown that kinematical constraints on the evolu- 
tion of non-dissipative quantum systems give rise to uni- 
versal bounds on the optimization of observables |2j and 
that the practical issue of dynamical realizability of these 
bounds depends on the controllability of the system ||. 
Controllability is also important in quantum computa- 
tion as it is directly related to the question of universality 
of a quantum computation element Qj. 

It has been shown that an atomic system with N ac- 
cessible energy levels, which are sufficiently separated 
to allow control based on frequency discrimination, is 
completely controllable using a sequence of frequency- 
selective pulses that address only a single transition at a 
time ||. Although this approach is enormously useful, 
control based on frequency discrimination is not appli- 
cable to systems with equally spaced energy levels such 
as truncated harmonic oscillators, and problematic for 
systems with almost equally spaced energy levels such 
as Morse oscillators, since in this case any external field 
is likely to address multiple transitions simultaneously. 
Moreover, systems with degenerate energy levels also 
present problems for this technique. 

In this paper we therefore concentrate on complete 
controllability of quantum systems subject to a single 
control field, e.g., a laser pulse or magnetic field, that 
addresses multiple transitions concurrently. In particu- 
lar, we study the problem of controllability of a system 



subject to a single control, for which the interaction with 
the control field is determined by the dipole approxima- 
tion. The systems considered are assumed to be non- 
decomposable, i.e., systems that can be decomposed into 
non-interacting subsystems are excluded. Note that such 
systems can never be completely controllable ||. 

It will be shown in particular that a non-decomposable 
quantum system with dipole interaction is completely 
controllable with a single control pulse if there is some 
anharmonicity in the energy levels. It must be noted 
however that it is sufficient if, e.g., the transition fre- 
quency for the first transition is different from all the 
other transition frequencies. This is a much weaker con- 
dition than frequency discrimination, which requires that 
all the transition frequencies are sufficiently different. 
Furthermore, for a system with no anharmonicity, i.e., 
equally spaced energy levels, we demonstrate that com- 
plete controllability depends on the values of the tran- 
sition dipole moments of the system, and establish suf- 
ficient criteria for complete controllability. The control- 
lability of some systems with degenerate energy levels 
is also discussed and examples of systems that are not 
completely controllable are presented. 



II. QUANTUM CONTROL SYSTEM 

Given any TV-level quantum system, the Hamiltonian 
of the unperturbed system can be written as 
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where {\n} : n = 1, ■ ■ ■ , N} is a complete set of orthonor- 
mal eigenstates and E n are the corresponding energy lev- 
els of the system. 

The application of external control fields perturbs the 
system and gives rise to a new Hamiltonian H = Hq+Hi, 
where Hi is an interaction term. In the control-linear 
approximation the interaction term is of the form 
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where f m (t) for m = 1, • • • , M are independent control 
fields, and the operator H m represents the interaction of 
the field f m (t) with the system. The off-diagonal ele- 
ments of H m depend on the transition dipole moments 
d n ,n' for transitions between energy eigenstates. Each 
H m is Hcrmitian since the transition dipole moments sat- 
isfy d n ' tn = d* nn ,, where d* n , is the complex conjugate 
of d n ^ n i . In this paper we are particularly interested in 
the case M = 1, i.e., a single control pulse, for which the 
interaction operator is of dipole form 

6 1 =^d n (\n)(n+l\ + \n+l){n\), d n ^0. (3) 

71=1 

Note that we exclude systems for which any of the tran- 
sition dipole moments d n = d n> n+i = d n +i } n vanish since 
these systems can be decomposed into non-interacting 
subsystems. 

An arbitrary initial state of the system can be repre- 
sented by a density matrix po that evolves according to 
the dynamical law 

p(t) = U(t,t )p U(t,t o y, (4) 

where U(t,to) is the time-evolution operator, which sat- 
isfies the Schrodinger equation 

ih^-U(t, t ) = (H + H!)U(t, to) (5) 
at 

with initial condition C/ (t , t ) = 1. 

We say the system is initially in a pure state if Tr(pg) = 
1. In this case the initial state of the system can also be 
represented by a normalized wavefunction IV'o)) which 
is either an energy eigenstate \n) or a superposition of 
energy eigenstates 

TV 

\i>o) = ^2 c n \n), (6) 

71=1 

where the c„ are complex coefficients that satisfy the nor- 
malization condition ^ n c„c* = 1 ||. The time evolu- 
tion of a pure state represented by a wavefunction \ip(t)) 
is 

m)) =U(t,toMo) (7) 

where \ipo) = \ip(to)) and U(t,to) is the time-evolution 
operator as defined above. 



III. CRITERIA FOR COMPLETE 
CONTROLLABILITY 

Since Ho + Hj is Hermitian, (^) implies that the time 
evolution operator U(t,to) is unitary. Hence, examina- 
tion of (Q) reveals that only target states p(tp) that are 
related to the initial state po by p{tp) = U PoU\ where 



U is a unitary operator, are kinematically admissible. 
However, in general, not all of these states can actually 
be dynamically reached, unless the dynamical Lie group 
generated by iHo and iH m , m — 1, 2, • • • , M, is the uni- 
tary group U(N). (See appendix ^ for a discussion of 
this requirement.) This motivates the 

Definition 1 A quantum system H = Ho + Hi with Ho 
and Hj as in ^ and (|j) is completely controllable if 

every unitary operator U is accessible from the identity 
operator i via a path j(t) — U(t,to) that satisfies (j^j. 

Complete controllability implies that any kinemati- 
cally admissible target state can be dynamically reached 
from the initial state by driving the system with a suit- 
able control field. If the system is initially in a mixed 
state represented by a density matrix po then this means 
that any other kinematically admissible mixed state can 
be dynamically reached. Similarly, if the system is ini- 
tially in a pure state represented by a normalized wave- 
function \ipo) then complete controllability guarantees 
that every other pure state represented by a normalized 
wavefunction \ipi) can be dynamically reached from the 
initial state. In |3j it is furthermore shown that complete 
controllability implies dynamical realizability of the uni- 
versal kinematical bounds on the optimization of observ- 
ables for (non-dissipative) quantum systems. 

It is apparent that if the dimension of the Lie al- 
gebra L generated by the operators {H , ■ ■ ■ , Hm}, 
or more accurately, their skew-Hermitian counterparts 
{iHo, • • • , iHm] is N 2 then Lo is the Lie algebra of skew- 
Hermitian NxN matrices u(N). Ramakrishna et al have 
shown in |?J , using results by Jurdjevic and Sussmann [|| , 
that in this case the dynamical Lie group of the system 
is the unitary group U(N). Noting that the dimension of 
u(N) is iV 2 and that any Lie algebra of skew-Hermitian 
NxN matrices of dimension N 2 is (isomorphic to) u(N) 
we have therefore 

Theorem 1 (Ramakrishna et al) A necessary and 
sufficient condition for complete controllability of a quan- 
tum system H = Ho + Hj with Ho and Hj as in (j^) and 
(Qj is that the Lie algebra Lq has dimension N 2 . 

This theorem provides a condition for complete control- 
lability of a quantum system that can easily be verified 
by computing the Lie algebra generated by Ho, ■ ■ ■ , Hm 
and determining its dimension. 

A basic algorithm for constructing a basis for the Lie 
algebra Lo in terms of iterated commutators is presented 
in Table Q. It can be optimized to increase the speed of 
the computation and to improve the accuracy of the nu- 
merical results. In Table | we use H to denote a N x N 
matrix and using the fact that a matrix can also be inter- 
preted as a vector, H for the A^ 2 column vector obtained 
by concatenating the columns of H vertically. W is a 
A 2 x TV 2 matrix whose columns W : j represent the basis 
elements of Lo- Note that W : j is the j-th basis element 
interpreted as N x A matrix. 
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let W = H 

let r — rank(W) 

for m = 2, • • • , M + 1 do 

if rank([W,H m ]) > r then 

append W by column vector H m 
r = r + 1 
endif 
endfor 
let r a = 
let r n = rank(W) 
repeat 

for / = r a + 1, ■ ■ ■ , r n do 
for j = 1, — 1 do 

let ft = [w;,;, W^] 

if rank([W,h]) > r then 

append W by column vector h 
r = r + 1 
endif 
endfor 
endfor 
let r = r n 
let r n = rank(W) 
until r n = r or r — N 2 

TABLE I: Algorithm to compute the Lie algebra generated 
by a control system {Ho, • • • , Hm} 

To initialize W we start with W- t i = Ho and add H m 
for to = 1, • • • , M provided that the additional column 
increases the rank of W. This guarantees that W ini- 
tially consists of linearly independent generators of Lq. 
To construct a basis for Lq, we compute all possible com- 
mutators of the columns W : j of W, interpreted as N x N 
matrices. Whenever a commutator is linearly indepen- 
dent of the columns of W, we add the commutator as a 
new column to W. Note that if we add a new column 
to W then we also have to compute the commutator of 
the matrix represented by the new column with all the 
matrices represented by the old columns of W. Hence, 
we repeat computing the commutators of the basis ele- 
ments represented by W until no new columns have been 
added in the previous step or the rank of W reaches the 
maximum of A^ 2 . 



IV. CONTROLLABILITY CALCULATIONS 

We implemented the algorithm presented in the previ- 
ous section and computed the dimension of the dynam- 
ical Lie algebra for the system H = Ho + f(t)Hi with 
Hq and Hi as in ([!]) and (||) for various choices of the 
energy levels E n and the transition dipole moments d n . 
In particular, we studied the AT- level harmonic oscillator 
with energy levels 

E n =»-§ (8) 

and the AT- level Morse oscillator with energy levels 

£„ = (n-i)[l-±B(n-i)] (9) 





Dim. N 


2 


3 


4 


5 


6 


7 


8 


Morse 


d n = y/n 


4 


9 


16 


25 


36 


49 


64 




d n = 1 


4 


9 


16 


25 


36 


49 


64 




d n = 1, djv_i = 2 


4 


9 


16 


25 


36 


49 


64 


Harmonic 


dn = s/n 


4 


9 


16 


25 


36 


49 


64 




dn = 1 


4 


4 


11 


11 


22 


22 


37 




d n = 1, djv-i = 2 


4 


9 


16 


25 


36 


49 


64 



TABLE H: Dimensions of the Lie algebra 



where B is a (usually small) positive real number. In 
our numerical computations we used B — 0.0419, which 
corresponds to a Morse oscillator model of the molec- 
ular bond for hydrogen fluoride We computed the 
dimension of the Lie algebra Lo for systems with varying 
dimension N and for different choices of the transition 
dipole moments d n . The results of some of these compu- 
tations are presented in Table ||. 

For the Morse oscillator system we observe that the 
dimension of Lq is always A^ 2 , i.e., it is completely con- 
trollable independent of the choice of the d n (as long as 
the d n are non-zero). 

For the harmonic oscillator, however, the dimension of 
Lq depends on the choice of the transition dipole mo- 
ments. For the usual choice, d n = y/n, the dimension 
of Lq is A^ 2 , i.e., the system is completely controllable. 
However, if we chose all the d n to be equal, e.g., d n = 1 
for n = 1, - ■ ■ ,N — 1, then the dimension of the Lie al- 
gebra Lq is less than A^ 2 and the system is therefore not 
completely controllable for N > 2. It is also worth not- 
ing that a slight modification of the d n is sufficient in 
this case to restore complete controllability: if we choose 
d n = 1, 77 = 1, ■ • • , N — 2 and (ijv-i = 2 then the dimen- 
sion of Lq is again A 2 . 

The extensive data we gathered strongly suggested 
that any Morse oscillator system with non-zero transi- 
tion dipole moments, i.e., d„ ^ 0, n = 1, • • • , N — 1, is 
completely controllable for any N, while complete con- 
trollability for a harmonic oscillator seemed to depend on 
the values of the transition dipole moments d n . These 
observations prompted us to study the issue of control- 
lability systematically using Lie algebra techniques. 

V. RESULTS FROM LIE ALGEBRA THEORY 

In order to prove our conjectures about complete con- 
trollability based on numerical evidence, a few general 
results from the theory of Lie algebras are required. For 
more detailed information about Lie algebras and Lie 
groups the reader is referred to @ 0, g |l| or any 
other book on the subject. 

We first observe that u(N) — su(N) "(1), where 
su(N) is the Lie algebra of traceless skew-Hermitian 
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N x N matrices. If the diagonal elements d m , m of the in- 
teraction operators H m for m > are zero, as is the case 
in the dipole approximation, then the interaction opera- 
tors H m are represented by traceless Hcrmitian matrices, 
i.e., iH m € su(N) for m > 0. If the internal Hamiltonian 
Hq is traceless as well, i.e., iH € su(N), then the dy- 
namical Lie algebra Lq generated by iHq and iH m must 
be su(N), or a subalgebra of su(N), since the commuta- 
tor of two traceless skew-Hermitian matrices is always a 
traceless skew-Hermitian matrix. By our strict definition 
of complete controllability, a system whose dynamical Lie 
algebra is su(N) is not completely controllable since its 
dynamical Lie group is SU(N), i.e., the Lie group of uni- 
tary N x N matrices with determinant one, and SU (N) 
is a proper subgroup of U(N), the Lie group of all unitary 
N x N matrices. (For a discussion of the practical sig- 
nificance of the difference between su(N) and u(N) see 
appendix |A|.) On the other hand, we have the following 
useful 

Lemma 1 If the dynamical Lie algebra Lq contains 
su(N) and Hq has non-zero trace then we have Lq — 
su(N)®u(l) ~ u(N). 

Proof: Note that we can write 

H = ±Tr(Ho)I + H , (10) 

where iH' € su(N). iH is in L by definition. Since 
Lq contains su(N) and hence iH , it must also contain 
the identity matrix il = (H - H' Q )N/ Tr(Ho)- Hence > 
noting that / generates a one-dimensional Lie algebra 
isomorphic to u(l) we have indeed Lq = su(N) © it(l). 
□ 

Thus, in order to show that a system is completely 
controllable we only need to show that Ti(Hq) ^ and 
that Lq contains su(N). 

To verify that Lq contains su(N) we need a complete 
set of generators for the Lie algebra su(N). Let e„,„' be 
the N x TV matrix such that the element in the n-th row 
and n'-th column is 1 while all other elements are 0, i.e., 

\&'n 1 n l )ij — 3inO~jn'> (H) 

where Sij is the Kronecker symbol. One can easily see 
that any traceless skew-Hermitian matrix must be a real 
linear combination of the N 2 — 1 basic matrices 

^n,n' = £n,n> - e n ',„ 1 < n < N — 1, n < n' < N 
^n,n' = K^n,n> + e n ',n) I < n < N - 1, n < n' < N 
h n = e nn - e n+ x, n +i 1 < n < N - 1. 

(12) 

However, verifying that Lq contains all of the N 2 — 1 
basis elements would be quite tedious. Fortunately, this 
is not necessary. 

Lemma 2 The skew-Hermitian N x N matrices e„n+i 
and e l n „ +1 , 1 < n < N, generate the Lie algebra su(N). 



Proof: Using the relation 

Cn,n'Cm,m' — Cn.m' '&n' 'mi 

which follows from the definition of e n ,ri > it can be veri- 
fied by direct computation that the skew-Hermitian ma- 
trices e^ n , and e ! n n , satisfy the equations 

«" + », (13) 

and a bit of algebra therefore shows that the 2(N — 1) 
elements e^ n+1 and e} n n+1 for 1 < n < N — 1 generate 
the entire Lie algebra su(N). □ 

Hence, if Lq contains the generators e^ n+1 and n+1 
for 1 < n < N — 1, then Lq must contain su(N). Thus, 
given a system whose energy levels are well enough sep- 
arated to permit selective control of each transition be- 
tween adjacent energy levels through frequency discrim- 
ination, i.e., 

N-l 

iH = iH + ^ fn{t) cos{^ n t)iHi 

n=l 

with H as in (0), /i n = E n - E n+1 and 

iH n = id n (\n)(n + 1| + \n + l)(n|) = d„e^ n+1 , 

where f n (t) is slowly time- varying compared to cos(/i„i) 
and d n ^ 0, we can conclude immediately that the Lie 
algebra of the system contains su{N) and hence that the 
system is completely controllable (if Tr(//o) 7^ 0). To 
see this, note that the generators e l n n+1 are given and 
the generators e^„ +1 can be obtained by computing the 
commutators 

V^-Q,e nn+1 ] = e^ Ln+1 . 

However, selective control of individual transitions be- 
tween adjacent energy levels through frequency discrimi- 
nation is not always possible. For instance, as pointed out 
earlier, it fails when the energy levels are equally spaced 
or degenerate, and it may not be a good approach for 
systems with nearly equally spaced energy levels, such as 
Morse oscillators. Furthermore, even if it is possible to 
use multiple pulses to selectively control individual tran- 
sitions, one may not wish to do so. Instead, one may 
for instance wish to control the system with a single op- 
timally shaped control pulse obtained using an efficient 
optimal control algorithm |lj, |lj, [l(| 

In order to establish criteria for complete controllabil- 
ity of iV-level systems subject to a single control field that 
drives all permitted transitions concurrently, we need an- 
other lemma, which makes use of the dipole form of Hi . 

Lemma 3 If H\ has the special form J9J, i.e., 

N-l 

iHt = dnei.n+i (14) 
n=l 
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then it suffices to show that Lq contains the pair of gen- 
erators ef 2 and e[ 2 or e^_ 1 N and e I N _ 1N , i.e., if Lq 
contains either of these two pairs of generators, then it 
contains all the generators of su(N). 

PROOF: If ef 2 ,e{ 2 G Lq then 

hi = 2 [^2^21] = *(en - e 2 a) G Lq, 

N-l 

Vt = iH x - d x e[ 2 = i^, n+1 G Lq. 

n=2 

where e 21 = e{ 2 . This leads to 

[hi,Vi] = d 2 e^ G Lq, 
-[hi, [hi,V{\] = d 2 ei 3 G Lq. 

Since d 2 7^ by hypothesis it follows that e 23 ,e 23 G Lq. 
Repeating this procedure N — 2 times shows that all the 
generators e^„ +1 and n+1 for 1 < n < N are in Lq. 
Similarly, we can show that e^_ 1 N and e I N _ 1 N in Lq 
implies that Lq contains all the generators e^ n+1 and 
ei +hn , l<n<N. □ 

Thus, if Hi has the special form (||) then it suffices to 
show that e.121 &12 € Lq in order to conclude that Lq is at 
least su(N). 



VI. COMPLETE CONTROLLABILITY FOR 
ANHARMONIC SYSTEMS 

Let fj, n = E n —E n+ \ for n = 1, • • • , N — 1 and V = iH%. 

Theorem 2 // /ii 7^ and fj,„ 7^ /if /or n > 1 £/ien i/ie 
dynamical Lie group of the system H — Hq + f(t)H\ with 
Hq and Hi as defined in (Q) and is at least SU(N). 
If in addition Ty(Hq) =fi then the dynamical Lie group 
is U(N), i.e., the system is completely controllable. 

Proof: We evaluate 

JV-l 

[iH a ,V] = -J2^nd n e^ n+1 = -V', 

n=l 
N-l 

{iHq,V'] = Y M^ne^,„+i = V". 

n=l 

Using V and V" we obtain Vi E Lq, where 
Vi = V"-/j 2 n _ 1 V 

N-2 



After N — 2 iterations, we have Vn-2 G Lq where 



Repeating the previous steps for i£/o and Vi leads to 
V 2 <E Lq, where 

U 2 = -[iHodiHoA]]- I&-2V1 

n=l 



AT-2 



7V-1 



n=2 



Since by hypothesis rfiII^=2 1 (A i i ~ A*n) 7^ this means 
Vjv_2 = ef 2 G io an d noting that /xi 7^ we have 



e 12- 



V, 



11 

N-2 



. — [iH Q ,V N _ 2 ] 



Mi 



G L . 



The conclusion now follows from lemmas [j] and ||. □ 

This theorem shows that for anharmonic systems com- 
plete controllability does not depend on the values of the 
transition dipole moments d n (as long as they are non- 
zero) and in particular we have the following 

Corollary 1 A quantum system H — Hq + f{t)H\ with 
Hq and Hi as in ^) and (Qj, respectively, and E n as 
in ^), i.e., a Morse oscillator, is completely controllable 
for arbitrary non-zero values of the transition dipole mo- 
ments d n . 

It is worth noting that Theorem [2] also applies to some 
degenerate quantum systems. 

Example 1 The system H — Hq + f(t)Hi with 

Hq = diag{Ei,E 2 , • • • , E 2 }, E x £ E 2 . 

and Hi as in (Qj is completely controllable by theorem 
U despite the fact that energy level E 2 has multiplicity 
N — 1 . In fact, the proof of controllability is even simpler 
for this system as we have V" = n\die{ 2 after only one 
step. 

This example begs the question whether the system is 
still completely controllable if we choose 

Hq = diag{£x ,---,Ei,E 2 }, Ei^E 2 

N-l 

instead. It is obvious that this system does not satisfy 
the technical condition [ii 7^ 0. However, it can easily 
be shown that this system is completely controllable by 
modifying the proof. In fact, this example is just a special 
case of the following 

Theorem 3 If fi^-i 7^ and /i^ 7^ t l ' N _ 1 for n < N — 1 

then the dynamical Lie group of the system H = Hq + 
f(t)Hi with Hq and Hi as defined in ^) and ^), re- 
spectively, is at least SU(N). If in addition Tr(Ho) =/= 
then the dynamical Lie group is U(N), i.e., the system is 
completely controllable. 

The proof of this theorem is analogous to the proof of 
Theorem 0. 
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VII. COMPLETE CONTROLLABILITY FOR 
HARMONIC SYSTEMS 

The condition fi\ ^ //^ in the previous theorems ex- 
cludes any system with equally spaced energy levels, such 
as a harmonic oscillator, for which we have /zi = //2 = 
■■• = hn~i- We shall assume that [i\ ^ since for 
/Hi = the system is completely degenerate with only 
one energy level. 

In order to state our theorem, we need to introduce 
some technical parameters 

2d\-d\, n = l; 

<•„ = i 2dl - dl_ x - dl +1 , n = 2,- •• ,N- 2; (15) 
2d 2 N _ 1 - d%_ 2 , n = N-l. 

Observe that these technical parameters depend on the 
values of the transition dipole moments d n . 

Theorem 4 The dynamical Lie group for a quantum 
system H = H + f \t)H\ with H and Hi as in (jl[) and 
^) and N equally spaced energy levels 

E n =E l + (n-l)ii U n = l,---,N, m ^ 0, 

is at least SU(N) if the the parameters v n satisfy one of 
the following conditions 

1. v n v N -i forl<n<N - 2; 

2. v n £ vi for 2 < 7i < N- 1. 

If in addition Tr(Ho) ^ then the dynamical Lie group 
is U{N), i.e., the system is completely controllable. 

Proof: Let V — iH\. In this case the element 

N-l 



are in Lq. Starting with V^ and , we have 



V = -^ 1 [iH ,V] = J2 d ^n, n +l 
n=l 

is in Lq and its sum and difference with V give rise to 



N-l 






n+1 + ^n.n+l) 


n=l * 


v ' 




N-l 






- e R \ 


n=l " 


v ' 



which, along with their commutator 
1, 



V? = ^i + ,*T 



N-l 



= id\ei.i + 2J i( d n - d n-i)^n,n - ■id 2 N _ 1 e NtN , 



N-l 



V+ = [V 1 °,V 1 + }-v 1 V+ 

N-l 



n=2 



JV-l 



= E( u ™ _U 2)K -vi)d n e+ n+1 e L , 

n— 3 



= d 



N-l 



N-2 
N-2 

Y[ ( V N-1 ~ V n ) 

71=1 



&N-1,N G L 0- 



Since by hypothesis Ujy— l ^ v n for n = 1, 2, • • • , TV — 2, 
we have e^-_ 1 jy € L . 

Similarly, starting with I/ and V{~ , we can also prove 
e JV _ 1 jy 6 Lq. This implies e^_ 1 w and e I N _ 1 N are in Lo 
and nence the conclusion follows from lemmas [l] and |^. 



71 = 2 



□ 



Example 2 The standard N-level harmonic oscillator 
system H = Hq + fit) Hi with Hq and H\ as in and 
E n = n — 1/2 and d n — yjn is completely controllable 
by theorem ^| since we have 

Vl = V-2 = ■ ■ ■ = V7V-2 = 0, VN-1 = N, 

i.e., the v n satisfy v^-i — v n = N for n — 1, ■ • • , N — 2. 



Example 3 However, the N-level harmonic oscillator 
system H = Hq + f(t)H\ with Hq and H\ as in and 
^j, E n = n — 1/2 and d n — 1 /or n = 1, • • • , N — 1 does 
not satisfy the hypothesis of theorem Q since we /iai>e 

= 1, V 2 = ■ ■ ■ = V N -2 = 0, V N -i = 1, 

i.e., vn-i = vi. Numerical calculations for various N 
confirm that this system is indeed not completely control- 
lable. (See Table 

For N — 4 one can, e.g., verify either numerically 
or by analyzing the Lie algebra that the simple skew- 
Hermitian matrix 



I % o o\ 

s/ i 

12 ~ 

V o o o o/ 
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is not in Lq and therefore the unitary operator 



U(9) = exp(0ef 2 ) = 



/ cos(0) isin(0) o\ 

ism(6) cos(0) 

10 

\ 1 / 



can not be dynamically generated for this system for any 
8 S (0,27r). Thus, if the system is initially in state po — 
Yln=i w n\ n ) ( n \ then the target state 

p(t F ) = u(d)p u(ey 

is generally not dynamically accessible since it is impos- 
sible to put the energy eigenstates |1) and |2) into super- 
position without equally "entangling" the states |3) and 
1 4) of the initial ensemble. 



Example 4 The N-level harmonic oscillator system 
H = H Q + f(t)Hi with H and Hi as in ^) and ||), 
E n = n-1/2, d n = lforn = l,-- -,N-2 butd N -i = 2, 
however, does satisfy the hypothesis of theorem Q since 
we have 



V\ = 1, v 2 



: VN~3 = 0, V N . 



-3, v N . 



7, 



i.e., v n 7^ v\ for n ^ 1. Numerical calculations for vari- 
ous N confirm that this system is completely controllable 
(see Table^). 

The rather surprising results of the previous two ex- 
amples can be understood by analyzing the Lie algebra 
generated by Ho and Hi in both cases. Although the 
details of this analysis are beyond the scope of this paper 
(and will be discussed in a future paper) we would like 
to mention here that for a harmonic system with transi- 
tion dipole moments satisfying the "symmetry" relation 
d-N-n = ±<^n for 1 < n < N/2, the transitions n — > n + 1 
and N — n—>N+l — n for 1 < n < N/2 become coupled. 
This leads to a collapse of the Lie algebra and loss of com- 
plete controllability, which can be restored by breaking 
the symmetry in the transition dipole moments. This is 
why changing d/v-i = 1 to d^^i = 2 restored controlla- 
bility in the last example. In fact, changing d^-i to any 
value other than ±1 would work as well. 



VIII. CONCLUSION 

The question of complete controllability of quantum 
systems using external control fields has been addressed 
before by various authors and it is, for instance, well 
known that a quantum system is completely controllable 
if it is possible to address a sufficiently large set of sin- 
gle transitions using multiple frequency-selective control 
pulses. However, many optimization strategies attempt 
to find a single control pulse that addresses all transitions 



concurrently to achieve the control objective. Further- 
more, control based on frequency discrimination is not 
always possible, e.g., it is not suitable for systems with 
equally or almost equally spaced or degenerate energy 
levels. 

Despite the relative importance of control strategies in- 
volving a single control pulse, sufficient criteria for com- 
plete controllability in this case have so far been missing. 
In this paper we addressed this problem and established 
general criteria for complete controllability of quantum 
systems subject to a single control pulse. 

In particular, we showed that most anharmonic, non- 
decomposable quantum systems are completely control- 
lable using a single control that drives all the transitions 
concurrently, independent of the values of the transition 
dipole moments d n . For quantum systems with equally 
spaced energy levels we demonstrated that complete con- 
trollability depends on the values of the transition dipole 
moments d n and derived conditions that guarantee com- 
plete controllability. 

We verified that the standard truncated harmonic os- 
cillator with transition dipole moments d n — i/n satisfies 
these conditions and gave examples of harmonic systems 
that do not satisfy the conditions. In the latter case 
we also checked by direct computation of the Lie alge- 
bra that they are not completely controllable and showed 
that there are certain unitary operators that can not be 
dynamically realized for these systems. 
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APPENDIX A: CONTROLLABILITY AND THE 
QUESTION OF TRACE 

All the theorems about complete controllability pre- 
sented in this paper require that the trace of H be non- 
zero. The mathematical necessity of this hypothesis is 
obvious: a set of traceless skew-Hermitian NxN matrices 
cannot generate all of u(N) but at most su(N). Hence, 
assuming that the interaction terms H m , m > 0, are rep- 
resented by traceless skew-Hermitian matrices, if Hq is 
traceless as well then the dynamical Lie group can be at 
most SU(N), the set of unitary matrices with determi- 
nant one, which is a proper subgroup of U(N). However, 
our definition of complete controllability requires that all 
unitary matrices be dynamically accessible. 

Nevertheless, the trace condition is physically some- 
what disturbing since the energy levels of a physical sys- 
tem are generally only determined up to a constant, and 
hence the trace of H seems physically rather insignifi- 
cant as one can always make it either zero or non-zero by 
shifting the energy levels of the system by a constant. We 
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shall attempt to resolve this apparent conflict by showing 
that the difference between SU(N) and U(N) is only a 
phase factor. 

Let the initial state of the system be represented by 
the normalized wavefunction |^o)' If the dynamical Lie 
group of the system is U(N) then any other pure state 
represented by normalized wavefunction \ipi) is dynami- 
cally reachable since given any two normalized wavefunc- 
tions there always exists a (not necessarily unique) uni- 
tary transformation U such that \ip%) = U\ipo) and we can 
find a path j(t) = U(t,t ) in U(N) such that j(t ) = 1 
and 7 (if) = U. Since the determinant of a unitary op- 
erator is a complex number of modulus 1, we can write 
det(7(») = e^W. Noting that det(aA) = a N dct(A) 
where N is the size of A, we see immediately that 
j(t) = e-^^jit) has determinant 1 

det(7(t)) = (e-^ (t)/JV ) Ar det( 7 (t)) = e-^We^W = 1, 
and thus defines a path in SU(N). Furthermore, 

7 (t F )|Vo) = e-**' N U\iM = e-^IVi), 

i.e., j(tF)\ipo) and differ only by a phase factor. 
Hence, if the dynamical Lie group is SU(N) then we 
loose control over the phase of the state, otherwise there 
is no difference. 

Thus, for practical applications that do not require 
phase control one need not worry about the trace. For 
instance, if the goal of controlling the system is to maxi- 
mize the expectation value of an observable A at a target 



time t F , 

(A(t F )) = m F )\Am F )), 

then clearly the phases of the target states are irrelevant 
as they are cancelled out by computing the expectation 
value anyway. 

Moreover, if the initial state is given by a density ma- 
trix po then any target state p(t F ) that is dynamically 
accessible via a path in U(N) is also dynamically acces- 
sible via a path in SU (JV). To see this, let j(t) — U(t, to) 
be a path in U (N) such that 

p(t F ) = U(t F , t Q )p Q U{t F , to)*. (Al) 

Again, we have det( 7 (i)) = e^ (t) and 7(f) — 
e"^ ( * )/JV 7 (i) defines a path in SU{N) that is equivalent 
to 7(4) since @ 

p(t) = 7 W/5o 7 (<) t 

= J(t)M(t) j , 

i.e., the phase factors e^ l ^^' N cancel out completely. 

However, there are some applications of control in 
quantum computation where it is important to have 
phase control and when SU(N) is not adequate. There- 
fore, we have chosen to require the dynamical Lie group 
to be U(N) for complete controllability. 
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